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Abstract 

At the electron-proton collider HERA the inclusive D*^ meson photoproduc- 
tion cross section has been measured with the HI detector in two different, but 
partly overlapping, kinematical regions. For the first, where {W^p) ~ 200 GeV and 
Q2 < 0.01 GeV^, the result is a{jp ccX) = (13.2 ± 2.2^?;^ The sec- 

ond measurement for < 4GeV^ yields cj(7p ccX) = (9.3 ± 2.1;;1:|| llg) A*^ ^'^ 
{W^p) « 142 GeV and a{-fp ccX) = (20.6 ± 5.51^;^ fib at {W^p) ^ 230 GeV, 
respectively. The third error accounts for an additional uncertainty due to the 
proton and photon parton density parametrizations. Differential cross sections are 
presented as a function of the D*^ transverse momentum and rapidity. The results 
compare reasonably well with next-to-leading order QCD calculations. Evidence for 
diffractive photoproduction of charm quarks is presented. 
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1 Introduction 



The study of heavy quark production in lepton-proton scattering provides an important 
testing ground for the standard model [Q. At the electron-proton collider HERA, heavy 
quarks are produced, according to QCD, by direct and hadronic (resolved) photon pro- 
cesses. The direct photon gluon fusion process 7(7 — > cc, where a photon emitted by the 
electron and a gluon from the proton generate a cc pair is expected to dominate. The 
major contribution is due to the exchange of an almost real photon (photoproduction), 
where the negative squared four- momentum transfer carried by the photon is ~ 0. 
The scattered electron is either lost in the beampipe or detected at small angles with 
respect to the electron beam direction. The fraction of cc events where the scattered 
electron is seen in the main detector (Deep Inelastic Scattering, DIS, > 4 GeV^j is at 
least one order of magnitude smaller 0. Measurements at HERA can be considered as a 
continuation of fixed-target photoproduction experiments 0, but at about one order of 
magnitude higher centre-of-mass (CM) energies, W-yp ~ C(200) GeV. 

In the Weizsdcker- Williams Approximation (WWA) [|| , the electroproduction cross 
section agp is expressed as a convolution of the flux of photons emitted by the electron, 
f^/e, with the photoproduction cross section 

aep = a{ep eccX) = J dy f^/dy) ■ cr{jp ccX), (1) 

where y represents the fraction of the electron energy transferred to the photon in the 
proton rest frame. For the direct photoproduction process, the cross section a^p, in turn, 
is assumed in leading order to factorize into the photon gluon fusion cross section and the 
gluon density in the proton 

a^p = a{'yp ccX) = J dxg Xgg{xg, f/) ■ a{'yg cc). (2) 

Here Xg denotes the momentum fraction of the proton carried by the gluon and fi the 
factorization scale. Estimates of the cross sections depend on the behaviour of the gluon 
density distribution g of the proton at small Xg, on the QCD renormalization scale, on 
the factorization scale, and on the heavy quark mass rric 0. 

Next-to-leading order (NLO) corrections to the parton cross section ^ are found 
to be substantial, but are reduced by experimental selection criteria, which limit the 
acceptance to the central rapidity range in the 7p CM system at large transverse mo- 
menta 0. 

Charm photoproduction can also proceed via the hadronic component of the photon 
(resolved photon processes), where a parton inside the photon scatters off a parton inside 
the proton, e.g. gg —>■ cc. This process known to dominate light quark production is 
expected to contribute much less to the production of charmed quark pairs |^ . However, 
the production cross section still depends strongly on the parton density function of the 
photon 0]. Other mechanisms, as for example the production of charm in the fragmenta- 
tion process, which is suppressed by the mass of the charm quark, or the production from 
the intrinsic charm content of the nucleon are believed to be small. These processes, 
as well as any possible intrinsic charm component of the photon, are neglected in the 
present analysis. 

Heavy quark production offers in principle the possibility of probing the gluon dis- 
tribution in the proton and the photon either indirectly, by measurement of the total 
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photoproduction cross section or of differential distributions, or directly, by the explicit 
reconstruction of Xg. Measurements of the first kind are described here. A similar mea- 
surement of (T^p has also been reported by the ZEUS Collaboration 

The analysis makes use of the D*-tagging [0], i. e. of the tight kinematical conditions 
in the decayQ D*^ D^n~^, where the mesons are reconstructed in the decay channel 
K^TT^. A better resolution is achieved in the distribution of the mass difference 

AM = M(L)°7r+) - M(D°) (3) 

than in the D*~^ mass distribution itself, whose width is dominated by the momentum 
resolution of the detector. 

The contribution of D*^ mesons, originating from decays of b flavoured hadrons is ne- 
glected, because of the expected small b production cross section at HERA (0(5 nb) ||TT[1 ). 
D*~^ mesons from decays of higher-mass charm states {e.g. D**) are not separated. 

Recently, much interest has been focused on a subclass of electroproduction events in 
which there is no hadronic energy flow in an interval of pseudorapidity, rj = — lntan(6'/2), 
adjacent to the proton beam direction. These diffractive processes are interpreted as an 
exchange of a colour-less object with the quantum numbers of the vacuum. The study of 
charm production in these processes is expected to provide information on the partonic 
structure of diffractive exchange. 



2 Analysis 

The present analysis is based on data collected with the HI detector during the 1994 
running period of the HERA collider, where 27.5 GeV positrons collided with 820 GeV 
protons, at a CM energy of 300 GeV. A detailed description of the detector and its trigger 
capabilities can be found elsewhere . 



2.1 Detector Description 

Charged particles are measured by two cylindrical jet drift chambers ||TB|, mounted 



concentrically around the beamline inside a homogeneous magnetic field of 1.15 Tesla, 
yielding particle charge and momentum from the track curvature in the polar angular 
ranged of 20° to 160°. Two double layers of cylindrical multiwire proportional chambers 
(MWPC) with pad readout for triggering purposes are positioned inside and in be- 
tween the two drift chambers, respectively. The tracking detector is surrounded by a fine 
grained liquid argon calorimeter |jT6|, consisting of an electromagnetic section with lead 



absorbers and a hadronic section with steel absorbers. It covers polar angles between 
4° and 155°. The luminosity is determined from the rate of Bethe-Heitler ep ep^ 
bremsstrahlungs events. The luminosity system consists of an electron detector and a 
photon detector, located 33 m and 103 m from the interaction point in the positron beam 
direction, respectively. The electron detector is used to tag photoproduction events by 
detecting positrons scattered at small angles. A time-of-fiight system (TOP) is located in 
the backward direction at z ^ —2 m. 



^Henceforth, charge conjugate states are always imphcitly included. 

^Hl is using a right-handed coordinate system with the z axis pointing in the direction of the proton 
beam (forward), the x axis pointing towards the centre of the storage ring. The direction of the incoming 
positron beam is termed backward. The polar angle is measured with respect to the proton direction. 
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2.2 Data selection and D*^ reconstruction 

The analysis is carried out independently for the case where the scattered positron is 
detected in the electron tagger and for the case where it is not required to be seen. 
Henceforth, the respective data samples will be referred to as tagged and untagged sample. 
They correspond to integrated luminosities of (2.77 ± 0.04) pb~^ and (1.29 dz 0.02) pb~^, 
respectively. About 20 % of the reconstructed D* candidates in the tagged sample are 
also present in the untagged sample. 

Proton beam induced background is reduced by requiring the event vertex to lie 
within ±40 cm of the nominal interaction point along the beam direction. A further 
reduction is achieved by excluding events with energy flow only into the forward region 
of the detector. 

For each event all possible M{K~tt~^) mass combinations are calculated with tracks 
of transverse momenta pt > 0.5 GeV/c. No particle identification is applied; each particle 
is assumed to be a kaon or a pion in turn. Pairs with an invariant mass within ±80 MeV/ c^ 
of the nominal mass of 1.865 GeV/c^ are combined with an additional track with 
Pt > 0.15 GeV/c and a charge opposite to that of the kaon candidate. 

Figure |1| shows the distribution of the mass difference (0) for D*"'"-candidates with 
Pt{D*+) > 2.5GeV/c and a rapidity -1.5 < y{D*+) = -|ln|^ < 1 for the tagged 
and untagged samples combined. D*^ production is found as a distinct enhancement, 
containing about 190 combinations in a ±2.5MeV/c^ window around the expected mass 
difference of 145.4 MeV/c^. No enhancement is observed if the mass difference for the 
wrong charge combinations M{K~7r~7r~^) — M{K~tt~) is used instead, as shown by the 
shaded histogram in figure |I[ The number of D* candidates is obtained from a simulta- 
neous fit to signal and background events in the right-sign (RS) and wrong-sign (WS) 
distributions of the AM spectra. The signal is described by a Gaussian and the back- 
ground shape is parametrized by a function of the form Oj- (AM — mjr)'', (i=RS, WS). The 
position and the width of the signal are determined from a fit to a larger data sample using 
additional trigger conditions, and then kept fixed at those values (AMq = 145.4 MeV/c^, 
(j(AM) = l.llMeV/c^) for all subsequent fits. Uncertainties from variations of the fit 
procedure are accounted for in the systematic error. 

A Monte Carlo simulation is used to determine the efficiency for the reconstruction, 
the selection cuts, and the acceptance of the detector. Hard scattering events for direct 
and resolved photoproduction are generated in leading order with the PYTHIA 5.7 pro- 
gram [|T^. The generated events are fed into the HI detector simulation program, and 
are subjected to the same reconstruction and analysis chain as the real data. 

The tracking efficiencies have been examined in detail using data. The single track 
reconstruction efficiency etrack is obtained by scanning tracks of high pt cosmic muons, 
where the measured pt of the incident track segment is compared with that of the outgoing 
track segment. The prdependence of etrack is determined directly from the data by a 
novel method [1^ based on the decay property of the pseudoscalar meson, decaying 
isotropically in its rest frame. The efficiency is found to rise from to the maximum value 
between pt = 90MeV/c and pt = 120MeV/c, and to remain constant beyond that. The 
precision measured in these studies is quoted as the systematic error. For single tracks the 
uncertainty found is ±2% for the track reconstruction and ^'^3% for associating the track 
to the primary vertex. Combining them for the three tracks yields a systematic error of 

on the tracking efficiency. 
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2.3 Analysis of electron tagged data 



Tagged events are required to have a positron candidate with energy Ee' > 4 GeV in the 
electron tagger and to have less than 2 GeV energy deposited in the photon detector. In 
addition, at least one charged track candidate has to be detected by means of a MWPC 
trigger |jl5|, [l^ and a drift chamber track trigger , thus ensuring activity in the central 



detector. The trigger efficiency is determined from the data itself, using independent 
triggers. The analysis is restricted to the kinematical region 0.28 < y = 1 — Ef.i / E^, < 0.65 
and < lO^^GeV^, where the acceptance of the electron tagger is above 20% with 
an average value of about 60 %. As a consequence, the 7p CM energy range is limited 
to 159GeV< W^p < 242 GeV, with a mean of W^p ^ 200 GeV and an average (Q^) ^ 
10"'^ GeV^. The efficiency excluding the y dependent electron tagger acceptance is found 
to be (48 ±4)% and (58 ±5)% for direct and resolved processes, respectively. 

2.4 Analysis of untagged data 

The untagged sample covers the kinematical region 0.1 < y < 0.8 and < 4GeV^ at 
an average (Q^) ~ 0.2 GeV^. Contributions from DIS with > 4GeV^ are rejected by 
requiring that no scattered positron candidate with E^' > 10 GeV be measured in the 
main detector. The remaining background from DIS events is suppressed by requiring 
y < 0.8 and is estimated to be less than 1 %. Here y is calculated from all measured final 
state particles using the Jacquet-Blondel method pi . 



The events are triggered by a combination of signals from the central and rear parts 
of the detector. At least one MWPC track candidate is required to point backwards into 
the region 110° < 9 < 155°. The backward TOF system must positively identify the 
event as a genuine ep collision by registration of a particle within the proper interac- 
tion time window and within its angular acceptance of approximately 160° < 9 < 177°. 
The trigger efficiency for the central MWPC and drift chamber trigger components is 
determined by simulation to be (84 ± 4)%. For the backward part it is obtained from 
data, imposing the same selection criteria but using independent triggers based on local 
coincidences of MWPC tracks and low threshold (> 1.5 GeV) signals in the liquid- argon 



(LAr) calorimeter ||2^. Sufficient statistical precision to determine the efficiency in bins 
of pt and y is achieved by including the sideband region of the mass difference signal, 
0.15GeV/c^ < AM < 0.18GeV/c^, and the wrong sign combinations. To account for 
the different event topology of the combinatorial background the efficiency is determined 
and parametrized as a function of y and then folded with the y spectrum of simulated D* 
events. This yields efficiencies of the backward trigger component of (28 ±4)% for direct 
and (35 ±5)% for resolved production processes, respectively. 



3 Results 

The visible production cross section in ep collisions is calculated from the observed number 
of D*^ mesons, A^, in the kinematical ranges of pt{D*) > 2.5 GeV/c and rapidity —1.5 < 
y{D*) < 1 according to the formula 

A^ 

aD'± = criep ^ D*^X) = a{ep ^ D*+X) + a{ep ^ D*-X) = ^ „ , (4) 

L ■ B ■ e 
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where L denotes the integrated luminosity, e the total efficiency, and B = B{D*^ 
L)°7r+) ■ B{D^ K-TT+) = 0.0273 ± 0.0011 is the combined branching fractions 
of D*^ and mesons. For the analysis of the tagged sample, the acceptance of the 
electron tagger and its trigger efficiency are accounted for on an event-by-event basis. For 
the relative ratio of direct to resolved photoproduction processes the values predicted by 
the NLO QCD calculation are used [i.e. 79:21 for the full {y,Pt) range or 93:7 for the 
visible kinematical range). The charm quark mass is assumed to be rric = 1.5 GeV/c^, 
the ratios of the factorization scale for the photon, the factorization scale for the proton, 
and the renormalization scale are taken to be 2mc, 2mc, and m^, as recommended by the 
authors ISll . 



3.1 Visible cross section a^p 

In the tagged sample the fitted number of D* mesons, 119 ± 16, is corrected for the 
electron tagger acceptance p4|, yielding = 197 ± 28. For the kinematical region 
PtiD*) > 2.5GeV/c2, -1.5 < y{D*) < 1, < 0.01 GeV^, and 159 < W^p < 242 GeV 
the visible ep production cross section is determined to be 

a{ep D*^X) = (4.90 ± 0.70i[J:^^) nb (tagged), (5) 

where the errors refer to the statistical and the experimental systematic error (see below). 

In the case of untagged events, the fitted number of D* mesons is 97 ± 15 events 
and the average total efficiency for 95 < W^p < 268 GeV and < 4 GeV^ is found 
to be 0.14 ± 0.03. The ep cross section in the kinematical region pt{D*) > 2.5GeV/c^, 
—1.5 < y{D*) < 1 is thus measured as 

a{ep D*^X) = (20.2 ± 3.3l|[!) nb (untagged). (6) 

The visible cross section is almost insensitive to both the choice of parton density 
parametrizations and to the mixture of direct and resolved photoproduction processes, 
because the efficiencies are very similar and there is no acceptance correction involved. 

The experimental systematic uncertainties are listed in table |1|. In the analysis of 
tagged data, the largest contribution (11 %) is due to the uncertainty in the track recon- 
struction, whereas in the untagged case the largest error arises from a 14 % uncertainty 
in the determination of the trigger efficiency. Adding the various uncertainties in quadra- 
ture results in a total experimental systematic error for the tagged sample of for the 
inclusive D* cross section and % for the charm cross section (see below), respectively. 
For the untagged sample the corresponding uncertainties are % and tH %. 

Predictions by the NLO QCD calculation for the visible cross section a£)*± in the 
tagged case assuming the following pairs of proton and photon parton densities of (GRV 
HO[0 + GRV-G HO[0), (MRSA'H + GRV-G HO), (MRSD0'|2^ + GRV-G HO) 



and (MRSA' + LACllH) yield values of 3.2, 2.8, 2.4 and 2.8 nb respectively. With the 
present accuracy a clear distinction is not possible, albeit a slightly better agreement is 
reached for parton densities with a rising gluon density distribution at low Xg. A similar 
conclusion results from the analysis of the untagged data. This is in agreement with 



measurements of the gluon density by other methods |^, e.g. from scaling violations of 
F2 in DIS. 
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Tagged 


Untagged 


Track trigger 


5 % 


5 % 


Electron tagger acceptance 


5 % 


— 


Backward trigger 




14 % 


Track reconstruction 


-Li 1 f\-f 


1 1 "I p,-/ 

% 


Signal extraction/background subtraction 


6 % 


6 % 


Luminosity 


1.5 % 


1.5 % 


D* , branching ratios 


4% 


4 % 


c ^ D* branching fraction 


7% 


7% 


Total experimental uncertainty 


-13/0 


-19/0 



Table 1: Experimental systematic uncertainties. 



3.2 Total cross sections and cr^^ 

The visible cross sections (within a hmited (y,pt) phase space) have to be extrapolated 
to the full {iiiPt) phase space to obtain the total cross sections. 

The individual acceptances for the direct and resolved processes as well as their 
relative strength depend on the choice of the parton densities, and therefore so also does 
the extrapolation performed by simulation. This is illustrated in table ^ which lists values 
for the acceptance calculated for various parton densities of the proton and photon and 
for different charm quark masses, for the kinematical region of the tagged sample. The 
numbers for the untagged case are similar. The acceptance is defined as the fraction 
of D* mesons within the quoted y and pt ranges with respect to the total number of 
produced D* . 

The derivation of the total cross sections is based on simulations using a charm 
quark mass of 1.5GeV/c^ and assuming the GRV LO parametrizations for both the 
proton and photon parton densities, which are in good agreement with measured parton 
densities. This leads to a charm ep production cross section of aep = (941 ± 160l|2o) 
at a/s = 300 GeV and < 0.01 GeV^ for the full y-range. The effect of hadronization 
is included using the fragmentation fraction Bc-^o*+ = 0.260 ± 0.021 [^. The ep-cross 
section is converted into a 7p cross section using equation (|l|), which yields for < 
0.01 GeV^ 

a(7p ^ ccX) = (13.2 ± 2.2tlj tH) jjh at {W^p) ^ 200 GeV. (7) 

For the untagged case the result over the range of 95 GeV < W^p < 268 GeV becomes 

a(7p ccX) = (12.6 ± 2.1+^-6 +9-4) at (W^p) ^ 180 GeV. (8) 

The third error indicates the additional extrapolation uncertainty as discussed be- 
low. The larger available kinematic range in W^p allows a division into the two regions 
95 GeV < W^p < 190 GeV and 190 GeV < W^p < 268 GeV, thus providing information 
about the energy dependence of the cross section. 

The results are summarized in table ^ for both analyses and compared in figure]^ 
with measurements by the ZEUS collaboration (at similar W^p) 0, and previous fixed- 
target experiments at lower energies 0. The inner error bars represent the statistical and 
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Proton 


Photon 


rric 






npirtoTi Hpimitv 


[GeV/c^l 


A pppnta n pp 


r^v>\r T r\ \tm\ 
LrKV LU [pUy 




1 O 

i.Z 


4.0 /o 


GRV LO 




1.5 


6.3 % 


GRV LO [0 




1.8 


10.8 % 


MRSG El 




1.5 


6.7 % 


MRSA' [H 




1.5 


10.4 % 


GRV LO [0 


GRV LO m 


1.5 


2.1 % 


GRV LO 


LACl |28 




1.5 


0.7 % 



Table 2: Acceptance for different parton density parametrizations for the direct or re- 
solved contributions, respectively, and for different charm quark masses, as used in the 
extrapolation from the visible to the total cross section. 



Quantity 


Tagged 


Untagged 


Untagged 


Range in W^p (GeV) 
Range in (GeV^) 


159 - 242 


95 - 190 


190 - 268 


< 10-2 


< 4 


< 4 


D* candidates 


119 ±16 


46 ±9 


51 ± 12 


{et,t)i%) (direct) 


29 ±2 


11 ± 1 


17 ±2 


Photon flux 


0.0141 


0.0486 


0.0155 


a{ep-^D*X) [y,Q'',pt,y] (nb) 


4.90 


11.0 


8.5 


Errors 


±0.70t°:I^ 






(7(7p ^ ccX) [y, Q2] (^5) 
Errors 


13.2 

_l_9 9+2.1 +9.9 


9.3 

_i_9 1 +1.9 +6.9 
=■=^•1-1.8 -3.2 


20.6 

p:+4.3 +15.4 
='=0.0_3 9 _7 2 



Table 3: Cross section results for tagged and untagged data samples. Errors shown are 
statistical, experimental systematic, and for a{'yp ccX) also uncertainties due to the 
dependence on the parton density parametrizations. 



experimental systematic errors added in quadrature. The outer set of error bars indicate 
the total error after adding in quadrature the extrapolation uncertainty discussed below. 
The cross section is rising by almost one order of magnitude as compared to the low 
energy measurements. 

Overlaid in figure ^ are predictions by the NLO QCD calculation with parton 
density parametrizations MRSG ||32| for the proton and GRV-G HO p5| for the photon. 
The upper and the lower solid lines delimit the range of values expected due to a variation 
of the renormalization scale within mc/2 < /i < 2mc. 

Calculating a^p with other combinations of parton density parametrizations in- 
creases the measured a^p by up to 75% (in the case of MRSG and LACl), or decreases 
a-yp by up to 35% (in the case of MRSA' and GRV-G HO). This variation reflects the 
uncertainties due to the choice of parton densities, and is quoted separately as a third er- 
ror in the cr^p cross sections. The extrapolation uncertainty due to fragmentation models 
(< 30%, estimated by a comparison with a cluster fragmentation as implemented in the 
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Herwig program [^]) and due to the choice of the charm quark mass (see table 2) is not 
included in the error. 

If the extrapolation is based on the proton parton density parametrizations MRSG, 
MRSA' or MRSDO' the value of a^p obtained (for the tagged sample at W^p = 200 
GeV) becomes (12.2 ± 2.0tli) /^b, (8.6 ± lAtH) fjh or (7.4 ± 1.2til) f^^b, which are to be 
compared with the QCD predictions of 9.8 fih, 6.0 /ib or 3.9 /ib, respectively. Measurement 
and prediction change in the same manner. Hence a total cross section measurement can 
presently not distinguish between the different gluon densities. 



3.3 Differential distributions 

Differential photoproduction cross section distributions, in the visible region where no 
model dependent uncertainties enter from extrapolation, are obtained by determining the 
acceptances and efficiencies bin-by-bin, separately for the two analyses. The distributions 
l/{2Bc^£)*+)-da{'yp D*^X)/dy are shown in figuresj^a and Ija integrated over the range 
2.5GeV/c < pt{D*) < lOGeV/c. The distributions 1/{2B^^d,+) ■ da{-fp D*^X)/dpt 
are presented in figures]^ and ^ for the rapidity range of — 1.5 < y < 1. The results from 
the analyses of the tagged and untagged samples are not combined because they cover 
different W^p (and Q"^) ranges. Note that the largest overlap between the two samples, 
namely 10 D* candidates, occurs in the bin of < y < 0.5 . The error bars represent the 
statistical error and, for the untagged data, the bin-by-bin systematic error due to the 
trigger efficiency. The other systematic errors of the overall normalization are identical 
to the errors quoted in table |I|. 

The histograms shown in figures ^ and ^ represent the absolute predictions of the 
QCD calculation including charm quark hadronization using the Peterson fragmen- 
tation function p4| (with parameter e = 0.06), and containing both direct and resolved 
photon processes. The calculations assume the parton densities MRSG for the proton 
and GRV-G HO for the photon, unless stated otherwise. The histograms are averages of 
calculations done at three representative W^p values, weighted by the photon flux inte- 
grated over the represented range. Good agreement within errors is observed between the 
shape of the measured pt distribution and the NLO QCD calculation. The y distributions 
on the other hand do not agree so well, with a possible excess of the data in the forward 
direction. 

To demonstrate the dependence on the charm quark mass, the predictions for masses 
1.2 GeV/c^ and 1.8 GeV/c^ are also given in figure |^ (dashed histogram). In figure ^ the 
influence of different proton parton density functions is illustrated by overlaying QCD 
predictions based on the MRSA' parametrization (dashed histogram). The effect of as- 
suming the LACl photon parton density parametrization (and MRSG for the proton) 
is marginal (dotted histogram). Although the total charm cross section is considerably 
larger when using LACl, most of the difference with respect to using for example the 
GRV-G HO parametrization lies in the forward region (2 < y < 4) and at low pt, outside 
of the visible range of this measurement. 
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4 DifFractive photoproduction of charm quarks 



A search for D*^ production by diffractive processes is performed in a sample of rapidity 
gap events, in which no final state hadronic energy flow is detected adjacent to the proton 
direction. The selection of diffractive events is based on a cut in rjmax < 2 and is described 
elsewhere [^] . Here rjmax denotes the pseudorapidity of the most forward calorimetric 
energy deposit in excess of 400 MeV. The selection of D* candidate events and the rejection 
of contributions from DIS are identical to those used in the analysis of untagged data as 
described above. However, because of the small signal expected, the analysis has not been 
restricted to a particular trigger condition, and thus comprises an integrated luminosity 
of 2.77pb-i. 

The mass difference distribution for the selected events is shown in figure exhibit- 
ing a clear D* signal. To substantiate the evidence for a diffractive production process, 
the rjmax distribution of the D* candidate events after background subtraction estimated 
from wrong charge combinations, is depicted in figure The distribution shows that 
most events have rjmax values close to the maximum possible value {rjmax = 3.65), but a 
clear excess of events with rjmax < 2 is observed. The data are compared with predictions 
of a non-diffractive model (as implemented in the PYTHIA program) and a hard 
diffractive model (RAPGAP ^6|). The former (latter) is normalized to the number of 
data events at rjmax > 3 {rjmax < 2). The sum of both models describes the shape of 
the data well (solid histogram in figure ^d), while the non-diffractive model alone fails to 
reproduce the shape of the rjmax distribution (hatched histogram). 

A lower limit on the visible diffractive cross section for the kinematical region < 
4GeV2, 0.1 < y < 0.8, pt{D*) > 2.5GeV/c, -1.5 < y < 1 and rjmax < 2 is derived. 
Assuming the trigger efficiency in this range to be 1 yields a conservative limit of 

a{ep D*^X) > 145 pb at 90% C.L. {rjmax < 2). (9) 



The cross section limit can be compared with the predictions of the diffractive model |^ 
which assumes a partonic structure of the diffractive exchange. If the diffractive exchange 
is dominated by a hard gluon |^ at an initial scale of Ql ~ 4GeV^, a cross section 
of 780 pb is predicted. On the other hand, if a quark dominated structure is assumed, 
the prediction is 29 pb in this model. The measured cross section is much higher then 
the latter prediction. Therefore, the data clearly disfavour a quark-dominated diffractive 
exchange. 



5 Conclusions 

Charm photoproduction cross sections have been measured through the detection of D*^ 
mesons. At an average 7p CM energy of 200 GeV the result is o"(7p — > ccX) = (13.2 ± 
2-2l?:7l4:8) /^b for Q2 < o.Ol GeV^ with (Q^) ~ GeV\ For the range < AGeV^ 
with '(g2) ~ 0.2GeV2 the values are a{-fp ccX) = (9.3 ± 2.ll||l|^)/xb at (W^p) ^ 
142 GeV, and (20.6 ± 5.51^^ 1^^2^) /ib at {W^p) ^ 230 GeV. These values are about one 
order of magnitude larger than those measured at previous fixed-target photoproduction 
experiments. Both the W^p dependence of the photoproduction cross section and its 
dependence on pt - and to a lesser extent on y - of the D* meson are reasonably well 
described by NLO QCD calculations. A slightly better agreement is reached with a steep 
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gluon momentum distribution in the proton. This is in accord with measurements of the 



gluon density by other methods pOf. The measured visible cross sections appear to be 
somewhat higher than the absolute QCD predictions. Evidence for charm production 
is found in a subsample of events which show a distinct gap of energy flow close to the 
direction of the proton and which can be interpreted as photon diffractive dissociation. A 
quark dominated diffractive exchange is clearly disfavoured by the present measurement. 
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0.14 0.15 0.16 0.17 

Am (GeV/c^) 

Figure 1: Distribution of the mass difference AM = M(ir~7r+7r+) — M{K~'k'^) for the 
combined tagged and untagged sample. The sohd dots represent the data, the hatched 
histogram indicates the background as obtained from wrong charge combinations. The 
sohd hue is a fit of a Gaussian plus a term for the background, as described in the text, 
with fixed width and position of the signal. 
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Figure 2: Total charm photoproduction cross section as a function of W^p. The sohd dots 
and stars represent the present analyses with statistical and systematic errors added in 
quadrature (inner error bars). The outer error bars indicate the total error if in addition 
the uncertainty due to the choice of parton density parametrizations is added in quadra- 
ture. The crosses refer to the results of the ZEUS collaboration, the other symbols indicate 
earlier measurements at fixed-target experiments. The solid lines represent the prediction 
of a NLO QCD calculation using the MRSG and GRV-G HO parametrizations of the 
proton and photon parton densities, respectively. The upper and lower lines delimit the 
range of values expected from varying the renormalization scale within 0.5 < ii/rric < 2. 
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Figure 3: Differential cross sections for the tagged sample (solid dots), (a) l/{2Bc^o*+) ■ 
da{-fp D*^X)/dy for events with pt{D*) > 2.5 GeV/c and (b) 1/(25^^^5*+) ■ da{-fp 
D*^X)/dpt for events with -1.5 < y{D*) < 1. The solid histogram shows the NLO QCD 
prediction, using the MRSG proton parton density parametrization with a charm quark 
mass of 1.5 GeV/c^. The upper (lower) dashed histogram indicates the effect of changing 
the charm quark mass to 1.2 (1.8) GeV/c^. The histograms are averages of calculations 
done at three representative Wyp values, weighted by the photon flux integrated over the 
represented range. Common systematic errors of (9(15%) are not shown. 
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Figure 4: Differential cross sections for the untagged sample (solid stars). (a) 
1/{2B^^D*+) ■ da{-fp D*^X)/dy for events with ptiD*) > 2.5 GeV/c and (b) 
1/{2Bc^j:i*+) ■ da{jp D*^X)/dpt for events with -1.5 < y{D*) < 1. The histograms 
show NLO QCD predictions for various parton density parametrizations for the proton and 
the photon: MRSG + GRV-G HO (sohd), MRSA' + GRV-G HO (dashed), and MRSG 
+ LACl (dotted). A charm quark mass of 1.5GcV/c^ is used for the calculations. The 
histograms are averages of calculations done at three representative W^p values, weighted 
by the photon flux integrated over the represented range. Common systematic errors of 



C(15 %) are not shown. 
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Figure 5: (a) Distribution of the mass difference AM = M{K~7r~^7r~^) — M{K~7r~^) for 
events with a rapidity gap with rimax < 2. The sohd dots represent the data, the hatched 
histogram indicates the background as obtained from wrong charge combinations. The 
solid line is a fit of a Gaussian function for the signal plus a background term, (b) i]max 
distribution of D* candidate events. The hatched histogram shows the prediction of a 
non-diffractive model (PYTHIA), the solid histogram the sum of the non-diffractive and 
a diffractive (RAPGAP) model. The former (latter) is normahzed to the number of data 
events at r]max > 3 {iimax < 2). 



